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Abstract 

We suggest that the new regime of J/ip suppression in Pb-Pb col- 
Usions found by the NA50 experiment at CERN is the result of non- 
trivial space-time evolution due to specific behavior of the Equation of 
State (EOS) near the QCD phase transition. We also study another 
suppression channel, the conversion of J /ip into r/c during the late cool 
hadronic stage, and find it rather inefficient. 
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1. Recent studies by the NA50 collaboration in Pb(158 GeV-A)-Pb 
collisions show significant J/V' suppression compared to an extrapolation of 
known trends in lighter beams and p-A collisions. The data show a new 
suppression mechanism for impact parameters b < 8/m. Surprisingly, this 
new suppression regime seems to set in very sharply. Such abrupt behavior 
may be a manifestation of dramatic collective effects, related to qualitative 
changes in hot/dense hadronic matter due to the QCD phase transition. Q 

Logically speaking, additional J/ ip suppression can be explained either by 
(i) increasing absorption rates or by (ii) increasing the time that J /ip spends 
in the hot/dense medium. In contrast to the available literature devoted to 
the first scenario we study the second one. An increase in the Quark-Gluon 
Plasma(QGP) lifetime was found in the hydrodynamic framework in [0, as 
a result of a discontinuity in the Equation of State (EOS) at the QCD phase 
transition. 

2. The oldest and the simplest argument for suppression is that once 
QGP is formed, free gluons can rather easily "ionize" J /if) in an elementary 
process gJ/ilJ —>■ cc similar to the photo-effect p, |^, |^. Later |^ it was 
argued that successive charmonium levels ip',Xj J/i^ "melt" at high temper- 
atures and densities due to Debye screening in the QGP, and ionizing gluons 
are superfluous. 

There is a signiflcant phenomenological difference between these two mech- 
anisms. "Melting" implies the existence of certain thresholds as a function 
of matter energy density and can therefore generate a rapid onset of suppres- 
sion. Furthermore, it naturally explains why different charmonium states are 
suppressed at different densities. So, although there are currently no realistic 
models of this type, this scenario may be successful. 

Gluonic excitation on the other hand, is simply a probabilistic process, 
and like any other absorption processes ("co-movers" for example) cannot 
create discontinuous behavior by itself. Therefore a recent paper proposed 
that the formation of QGP is discontinuous; QGP is produced suddenly after 
a certain critical energy density ec is reached. This work was motivated by 
the classical nucleation theory of "bubbles" at a flrst-order phase transition 
and has a number of theoretical and phenomenological problems. First, the 
thermodynamics of the transition tells us that plasma formation is discon- 

^ We stress however, that we discuss the new trend over the whole region and not just 
its onset. 
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tinuous in temperature but not in energy density. A mixed phase linearly 
interpolates between the energy densities of the two phases, e^m and emax- Q 
Second, this scenario has a phenomenological problem: a jump into the QGP 
phase at some ec changes the equation of state discontinuously and implies a 
jump in entropy (or multiplicity) ||21||. Preliminary data [|l| on < ^^^(-Et) > 
shows no discontinuity. 

3. As deduced from lattice simulations at finite temperature, the QCD 
phase transition leads to a very non-trivial EOS. The ratio of pressure to 
energy density, p/e has a deep minimum at the end of the phase transition, 
^max ~ 1 — 2GeV / fm? . For the resonance gas and the high density QGP 
this ratio is relatively large (.2 and 1/3 respectively) but is only about .05 
at this minimum. This minimum is referred to as the softest point of the 



EOS and is associated with a slow hydrodynamic expansion, [g, |TT| When 
the initial energy density was close to the softest point, the QGP was found 
to live 2-3 times longer than the more energetic and more central collisions 
at the SPS.@ 

In this letter we study the extent to which the non-trivial features of 
J/ip suppression can be explained by such an increase in the QGP lifetime. 
The hydro solution in was made only for azimuthally symmetric central 
collisions. Assuming that about half of the collision energy is contained in 
the fully stopped part of the matter, the above authors expect emax to be 
reached at Piab ~ AOGeVA for central AuAu collisions. (As this lies between 
the nominal AGS and SPS energies, they proposed making additional low 
energy runs at the SPS.) Of course, the initial energy density decreases with 
impact parameter and since all values of b are triggered on anyway, scanning b 
is easier than scanning the collision energy in the experiment. Examining the 
dependence of energy density on impact parameter, one can see that indeed 
an impact parameter b ~ 8/m in Pb-Pb collision roughly corresponds to the 
energy density of the observed change of trend in J/ip suppression. 

Before turning to specifics, we point out a qualitative difference between 
our scenario and the one based on "melting". In the latter case the various 
charmonium states are completely suppressed for all e > ec, while in our 



^The arguments of |^ rely on the non-equilibrium kinetics of bubble formation. Note 
however, that a very small surface tension found by lattice studies, about O.OIT^^ makes 
the formation of small bubbles quite probable . 

^ Roughly, since p is the moving force and e is the mass to be moved, acceleration is 
proportional to this ratio. 
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case the suppression is centered around e ~ tmax and the survival probabihty 
increases for e > Cmax- This difference is manifest in the specific predictions 
discussed below. 

Hydro evolution for non-central collisions is much more complicated be- 
cause of a directed flow (see e.g. |13[)- Fortunately, we avoid these effects 
by restricting the discussion to J/tp suppression. For the first few fermi/c 
expansion is predominantly longitudinal, and thus similar for central and 
non-central collisions. 

4. It is clear that significant changes in the QGP lifetime should affect 
J/ip suppression, even if absorption rates are unchanged. However, J/ijj 
may escape from the system, which obviously limits the sensitivity to QGP 
lifetime. 

In order to reproduce many known observables, such as the distribution 
in the transverse energy Et and its correlation with b, we have constructed a 
small event generator to model our scenario. For Pb(158-GeVA)-Pb collisions 
at a given impact parameter, we generate charmonium events which have 
survived nuclear absorption and are now in a hot de-confined medium. The 
survival probability for a given event is e"^**" where tp is the time spent 
transversing this medium and P is the ionization rate. P is chosen to match 
the observed 25% suppression immediately following the jump. Our goal 
is to provide a reasonable explanation for the region near the discontinuity 
without too much regard for the most central points. 

In our model J/ ip interacts with the hot medium only when the wounded 
nucleon density, 71^, is greater than 3.1 fm~^. This condition is similar to [O 



and determines the transverse size of the region. We estimate the longitudinal 
size from the Bjorken formula. The average energy density reached in a 
collision at impact parameter b is 

3 dE° 3q 

< e >=< n TT- >= n < ^wib) > 

AycT dA AycT ^ ^ 

where < > is the average wounded nucleon density, and q = A GeV is a 

calorimeter dependent constant as defined in [jl5[. Since the jump is observed 



for b ~ 8.5 fm and < riwip = 8.5) >~ 2.6 fm~ , we choose r = 1.9 fm to match 
the energy density of the jump to the softest point, emax = 1-6 GeV/iw? . The 
thickness is then 2cr, or the distance the colliding pancakes have receded 
during the equilibration time. 
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In this model J / ip suppression begins only for wounded nucleon densities 
that are 20% larger than the softest point average, i.e, significant suppres- 
sion begins only above the mixed phase, e^ax • For the lifetime of this hot 
region which destroys J/ip we take the forms shown in Fig. la and calculate 
the corresponding suppression patterns shown in Fig. lb . To calculate the 
time inside plasma we need to distribute these events in coordinate and mo- 
mentum space. The distribution in the transverse plane is given by Glauber 



theory, like in |T5| and the distribution in the longitudinal direction is pre- 
sumed uniform. We have calculated transverse momentum distribution of 
J/ijj within Glauber theory (as in |T^). The longitudinal momentum or 
distribution is known from p-A data and in the plasma frame follows the 
form ~ (1 - |x/ 1)3-9. ^ 



We can then determine the time tp that J/ip remains in plasma and the 
corresponding survival probability, e~^^p, and find the J/ip yield 

' ''''-mf{z)hAPT)f{xf)e-^^^ 



db^ 

Here s is the coordinate in the transverse plane and each / is a distribution 
functions described above. Such a b dependent cross section may be converted 
to an Et dependent cross section using the Ex—b assignment of |T5[. Dividing 



this cross section by the Glauber value we compare our survival probability 
with preliminary NA50 experimental results in Fig. lb . 

The usual Glauber-type J/ip absorption on nucleons describes well p- 
A , lighter ion, and peripheral Pb-Pb data and and therefore we plot only 
deviations from this trend in Fig. lb . The increased QGP lifetime leads 
to additional suppression for intermediate Et between 50 — 100 GeV. Since 
J/ijj may leak out of the plasma region the discontinuity due to the lifetime 
is softened and we do not reproduce the very sharp jump seen in the data. 
However, agreement with most of the points in the intermediate region is 
significantly improved. Furthermore, our curves reproduce the vanishingly 
small slope of the suppression in this region. 

5. The pT-dependence of the suppression may separate the different sup- 
pression mechanisms discussed above. We have therefore calculated < > 
vs. Et- As a function of b, < Pt > is given by 
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Converting the numerator and denominator to Et dependent cross sections 
we plot < > vs. Et in Fig. Ic . As a result of the usual initial state 
parton re-scattering, < > increases with the centrality of the collision. 
If J/ip "melt" in the central region then only peripheral J/ip survive, and 
< Pt > decreases for the most central collisions, as suggested by Kharzeev 
et al. In our scenario however, some J/ip may still "leak" even from the 
central regions, and this flattens out the < p^ > vs. Et dependence. This 
dependence was indeed observed by the NA50 collaboration 




b(fmj ° " '° " " E, (GeV) ° " " " " " (GcV) 



Figure 1: For the different lifetimes of the plasma region shown in (a) the 
corresponding suppression pattern is shown in (b). (c)The J/ip pt distribu- 
tion of this work (curve C in (a)) compared to [16]. 

6. This part is rather independent from the rest of the paper; Here we dis- 
cuss the conversion of J/tp into its spin partner r/c. This channel is one more 
potential "sink" for J/ip, which (to our knowledge) has not been investigated 
previously. 

These reactions may be important in the later stages of the collision. The 
"photo-effect" above requires ~ 1 GeV gluons to ionize J/ip and therefore 
can only proceed during the hot plasma stage of the collision. In contrast, 
since J/ip and rjc are nearly degenerate, spin exchange reactions require little 
input energy and may proceed during the cool hadronic stage, which lasts for 
10-20 fm/c in Pb-Pb collisions. If these reactions have significant rates, some 
fraction of J/^p can be converted into rjc. The inverse reaction is of course 
also possible, but the rather short natural lifetime r^^ ~ lb fm/c prevents 
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significant feedback. 

Tlie two "spin excliange" reactions we consider liere are ip + n ^ rjc + p 
and ip + p rjc + TT. Below we calculate these reaction rates in a hadronic 
gas at T=150MeV using a non-relativistic quark model. 

First we consider the reaction ip + ir rjc + p. In a non relativistic 
quark model this reaction proceeds via a spin exchange collision between the 
c quark in ip and the u quark in vr. We use a point like spin-spin interaction 
between the two particles in the center of mass frame of the pion and J/ip. 

Hint = C(fc ■ (fdS^^r^ - r^) (1) 



From D-D* meson splitting Q we take the value C = 4.4 (GeV)'"^. Using 
this interaction Hamiltonian and assuming identical spatial wave functions 
for 7r,p and ipyfjc, one obtains the following cross section 

= ;7^K^V^ ■ / ?^^(2-)'5'(p-) (2) 

The unpolarized cross section is found by summing the cross section over 
final spins and averaging over initial spins 



q2 



EpPp 



{{a. a)) (3) 



where {{a ■ a)) = ^ 51 K^V^ • 

^ ss' 

for J/tp + 71 ^ r]c + p we have {{a ■ a)) = 1 and for J/tp + p—^ric + TTwe 
have {{a ■ a)) = 1/3. The value of the cross section(3) is determined from 
the incident energy of the pion and the kinematics of the collision. For a 
stationary J/ip the threshold for this process occurs when E„ ~ 710 MeV. 
The cross section is about 1.2mb at 800 Me l^. 

Now we calculate the reaction rate in an ideal pion gas for the process 
J/ip + 7T Tif. + p. The J /ip is at rest in the pion gas. In this energy range the 

^The exact physical nature of the spin-spin forces between hght and heavy quarks is 
not important for our results. Note however that it may be either due to the one-gluon 
exchange (as assumed in the works mentioned above) or due to small size instantons, see 
S. Chernyshev, M.A. Nowak, I. Zahed, Phys. Rev. D53, 5176 (1996). 
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pion mass is neglected and Boltzmann statistics are reasonable. The reaction 
rate is given by (flux) x (cross section) 

W^^,^ = 1^ / e-r^v,Ela{E,) dE, (4) 

/TT J threshold 

where gfTr = 3 is the isospin of the pion. Using the cross section above (3) we 
have calculated the reaction rates numerically 

Wj/^+^^r,.+p ^ .29 MeV (5) 

A very similar calculation gives the reaction rate for J/ip + p rjc + tt. 
For example the unpolarized cross section is given by C'^/tt'^ \p.,^E^/{vp + 
Vj/ip)]cm{{c ■ cr)). For this process we find 

Wj/^+p_.r,c+n ^ -^0 Mev (6) 

It has been suggested that the mass of the p shifts downward in medium, 
and the CERES dilepton data seem to support this claim. For reference we 
calculate the above rates for nip — 385 MeV 

Wj/^+^^r,c+p ~ (7) 

Wj/^+p^r^^+^ ^ 1.1 MeV 

Even these rates are rather small, and we conclude that it is unlikely that 
this channel can kill more than few percent of J/V'- 

Finally, we consider experimental checks of rjc production. Its 77 branch- 
ing is too small to be useful. Among the known modes the most promising 
seems to be a 00 decay with branching ratio of .7%. As the 00 background 
is not large, identification/mass resolution in the KK channel is generally 
good, and the rjc peak is rather narrow, rjc identification may be feasible. 

7. Discussion and Outlook. Recent findings by the NA50 collaboration 
indicate that an additional J/i/j suppression mechanism sets in for b < 8fm. 
This suppression can be created by a change in the absorption of J/ip, or by 
a change in the space-time evolution, or both. "Melting" may increase the 
absorption of J/ip in plasma, while a change in the EOS may increase the 
lifetime of the plasma itself. 

Two qualitative differences separate these scenarios experimentally. First, 
"melting" is not directly related to the phase transition, but is a subtle phe- 
nomenon where basically % states no longer exist. The space-time evolution 
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on the other hand, is directly related to the EOS and hence to the tran- 
sition point itself. Experimentally therefore, a change of behavior at the 
corresponding impact parameter Q in the hadronic observables sensitive to 
the transverse expansion , e.g. the pt slopes (especially of heavy secondaries 
N,d), HBT transverse sizes etc., would signal a change in the space-time 
evolution of the system. 

A second difference, is that while "melting" is presumed to destroy com- 
pletely the various charmonium states, the absorption in our scenario comes 
from the usual probabilistic rates, and therefore J/ip may escape from the 
central plasma region. This produces a qualitative difference in the mean 
transverse momentum vs. centrality, as discussed above. Much more detailed 
studies of the pt dependence of the suppression are now underway and should 
clarify this issue. 

There remain some unresolved experimental issues with the preliminary 
data we use. More statistics and further analyses are needed to explain some 
fluctuations in the data points and most importantly to verify the shape of 
the transition to the new regime, which now appears as a discontinuous jump. 
Finally, it is important to test experimentally an assumption, common to all 
explanations given above, i.e. that the phenomenon is induced at some value 
of energy density. This can be done by lowering the collision energy of the 
Pb beam. 
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